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NOMENCLATURE 
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cross-section area 
load damping coefficient 
leakage coefficient 
force 

load spring constant 
valve pressure coefficient 
valve flow gain 
flow force spring rate 
mass 

load pressure on power piston 

load pressure on input piston 

supply pressure on power piston 

supply pressure on input servovalve 

flow rate 

time 

volume 

valve area gradient 
input piston displacement 
valve displacement 
power piston displacement 
effective bulk modulus 
damping factor 
natural frequency 
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SUBSCRIPTS 
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refers 


to 
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refers 


to 



input piston 
load 

power piston 
piston 

supply pressure 
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I. 



INTRODUCTION 



During recent years, naval applications of gas turbine engine 
systems have become more prominent. At first, the gas turbine engine 
was developed as an aircraft power source. Presently, it is entering 
service as a marine power source. Gas turbine engine systems have 
several advantages compared to other power systems. The main advan- 
tages are its relative compactness and the fact that the rotating out- 
put of the turbine is usually what is wanted to turn a wheel, propeller 
generator, or pump. Unfortunately, compared to other systems, there 
is a problem of matching the characteristics of the load and the engine 
In this respect, controllable reversible pitch (C.R.P.) propellers are 
useful. 

Controllable pitch propellers, unlike fixed pitch propellers, have 
the ability to reverse thrust through a reversal of pitch rather than 
a reversal of rotation. The controllable pitch propeller is much more 
complex than a solid propeller and contains many more parts [1] . How- 
ever, the operation of a C.R.P. system provides a greater degree of 
flexibility in the maneuvering and docking of ships. 

Since there is an increasing interest in the application of con- 
trollable reversible pitch propellers to naval vessels, the objectives 
of this study were as follows: 

1. to determine an overall analytical model of an operational 
C.R.P. system; 

2. to check the system response for design variations. 
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II. BACKGROUND 



A. BASIC CONTROL CONCEPTS 

In a control system in which the output appears as mechanical energy , 
the quantity to be controlled can be either load, torque, position, or 
velocity and can be controlled under either open- or closed-loop condi- 
tions. In an open-loop system [2] , the input is fed to the controller 
which then operates on the output. The input, being an independent vari- 
able, is incapable of correcting any errors between the desired and 
actual output. 

The accuracy of an open-loop system is, therefore, greatly affected 
by external disturbances and variations of system parameters. In a closed- 
loop system [2] , the output is fed back and compared with the desired 
value. The error between the desired and actual output is fed into the 
control elements which in turn operate on the output, tending to reduce 
or null the error. This is the concept of negative feedback. A block 
diagram of a closed-loop system is shown in Figure 1. 




Figure 1. Block Diagram of a Closed-Loop System. 
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In a position control system in which the manipulated variables are 
pressure and flow, the hydraulic quantities are converted into a dis- 
placement by action on the piston rod or output shaft. The two proper- 
ties of flow which can be manipulated are its direction and magnitude. 
Directional control of flow is used for simple position control systems 
under both open- and closed- loop conditions [3] . 




1 Ps 

Supply 



Return 



Figure 2. Valve Controlled Piston. 

Directional control valves that could be used here are of the three- 
or four-way type. A schematic arrangement of a four-way valve control- 
ling a symmetrical piston is shown in Figure 2. 

When the control valve as shown is in the neutral position, the 
piston is hydraulically locked. A displacement of the valve spool in 
either direction causes the piston to move in the opposite sense. 
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An idealized graph of a position control system which is operated 
under open-loop conditions, by a directional control valve is shown in 
Figure 3, 




Figure 3. Position Control System Response (Open- Loop) . 

When the valve is displaced from its neutral position, the actuator 
moves at constant velocity in the required direction until the valve 
is returned to neutral. An open- loop system of this type is only suit- 
able for simple applications which do not require accurate piston control. 

Closed-loop position control systems employ controllers having non- 
continuous or "on-off” characteristics. When the actuating signal 
reaches a threshold value, the directional control valve becomes opera- 
tive and causes the actuator to move at constant velocity in the chosen 
direction. It continues to move until the error signal falls below the 
threshold value* This causes the valve to return to its closed position 
and stop the motion of the piston. 
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The most widely used type of control for hydraulic positioning servo 
systems is the control flow type, where both flow direction and magnitude 
are manipulated. The two methods of achieving flow control are: 

1. controlling the displacement of a variable delivery pump; 

2. throttling fluid through the variable orifices of a control 
valve. 

The two basic differences between pump operated and valve operated 
servo systems can be described as follows: 

1. In pump operated servos, working pressure is variable, 
depending on the external loading, while valve controlled 
servos operate at constant system pressure. 

2. In valve controlled servos, power sources are outside 
the loop. 

In this study, a valve controlled servo system was studied. 

B. ANALYTICAL APPROACH 

This analytical investigation was directed toward predicting various 

C. R.P. system pressures and motions. An analytical model was developed 
using first-order differential equations. The mathematical model, which 
was programmed on a digital computer, predicted the position and velocity 
of the power piston and the pressure at various points in the system as 

a function of time. Further analysis was performed using the model to 
obtain preliminary predictions of system performance under a variety of 
external loadings . 

To understand the derivation in the following sections. Figure 4 
indicates the geometry of a basic controllable reversible pitch (C.R.P.) 
propeller system. Beginning with the input side, a force balance on the 
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input servovalve yielded an equation for the load pressure on the input 
piston in terms of the valve motion. The input piston was moved by load 
pressure, P n ^. The valve rod control assembly was moved by the force 
exerted on the input piston. This caused the primary supply pressure, P^, 
to pass through the hub valve into the power piston and move the propeller 
crank-slot pin. This power piston control system, as such, is similar in 
operation to a vehicular power-assisted steering mechanism. 

To facilitate the digital simulation it was assumed that the input 
valve spool displacement, X^, and the arbitrary load, F^, were step 
inputs. In the process of solution, the system state equations were 
solved using CSMP III, the IBM simulation language. 
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III. CONSTITUENT RELATIONSHIPS 



It was a difficult problem to identify the dynamic characteristics 
of the controllable pitch propeller system because of its essentially 
nonlinear behavior. Identification was partially achieved by assuming 
that these responses were linear and then applying parameter estimation 
for linear systems. 



A. INPUT SERVOVALVE AND PISTON 

The combination of a servovalve and a piston is a common hydraulic 
power element [3] . In analyzing this simple system, three basic rela- 
tionships surface. The first concerns the valve, the second concerns 
the piston dynamics and the third concerns the flow of hydraulic fluid 
between the elements. 

Regarding the valve, the basic equation that pertains is 



where 



Q 0 = K X - K P, 
ica qa va ca la 



K 

qa 



K 

ca 



70 w v4? 
a sa 




P 

sa 



( 1 ) 



Regarding the piston, the basic equation governing value dynamics is 

ZF = M X (2) 

P v 
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The sum of the force acting on the piston is 



ZF = F „ - F n - K X-B X 
p pic Jca pa pa 



( 3 ) 



where 



F = A P 0 

piC p Jca 



F n = K- (X - Y) = 0.86 w (P - P ) (X - Y) 
iCa f m s m 



Combining equations (2) and (3) gives 



A P 0 = M X + B X+K X+F, 
p ica v pa pa ica 



(4) 



The flow relationship between the two elements is developed by 
applying the continuity equation to each of the piston chambers. In 
final form, this continuity relation becomes 



V 

Q £a = C tp£a P £a + A p * + ( 4 6 ) 

e 



(5) 



where is an overall leakage coefficient. 

Equations (1) , (4) and (5) are the basic equations of this sub- 

system. 

B. POWER PISTON 

From the characteristic equation of the critical center spool valve 
(hub valve) , one can derive 



Q = K (X-Y)-K P 
m qm cm m 



(6) 



where 



K 



qm 



140 w /p 
m s 



K 



'em P 
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In a manner similar to the development of equation (5) , the continuity 



expressions for the hub valve can be combined to yield 



V 

• m 

0 = C . * P + A Y + (- 7 - 5 - 

~m tp£m mm 4 6 



-) P_ 



m 



(7) 



where C ^ is an overall leakage coefficient. 

The final equation arises by applying Newton’s second law to the 
power piston [4] . The resulting force equation is 



A 

m m 



m £ y 



+ B 



pm 



(X - Y) 



+ K (X - Y) + F p 
pm x, 



( 8 ) 



Equations ( 6 ) , (7) and ( 8 ) are the basic equations of the power 

piston system. 

C. OVERALL SYSTEM 

In analyzing the overall system, two modes of operation became 

apparent. First, if the system input was a forward valve displacement, 

X , the output was a power piston position and velocity. The second 

mode occurred during recovery from an external loading. For this mode, 

the input was an external load, F^. At this point, it was necessary 

to simplify the governing equations. Usually, the loading on a valve 

controlled piston system is simpler than the most general in two 

respects [3] . First, significant spring loads are usually absent and 

hence K = K = 0. Second, damping loads are likewise negligible and 

then B = B =0. 
pa pm 



20 



Thus, equation (8) reduced to 



Y = (A P - F„) 

mm % 



( 9 ) 



A combination of equations (6) and (7) produced 



4 6 e 

? = ^ (- K P-AY + K X-K Y) 

m V cem mm crm qm 

m 



( 10 ) 



Equation (4) was simplified to 



x * — <- K f X + K f V + A p P la> 

V 



( 11 ) 



Finally equations (1) and (5) yielded 



Za 



4 6 



V 



(- K 



cea 



Za 



-A X + K 



qa 



X ) 
va 



( 12 ) 



Equations (9-12) were the governing equations for the C.R.P. system. 
The following is a listing of the constants used in the derivation of 
these governing equations and in the accompanying computer solution. 

The valves were derived from manufacturer's data [5]. 
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TABLE I 



A (AA) 

P 

A (AM) 
m 

K (CEI1) 
cm 

K (CEA) 
ca 

V (VA) 
a 

V (VM) 
n 

w ( — ) 

cL 

w (W) 
m 

P (PS) 
s 

P (PSA) 
sa 

S e (BE) 

( XM) 



M v (VM) 



K (QAK) 
qa * 



K 

qm 



(QMK) 



LIST OF CONSTANTS 



cross-sectional area of input piston 
cross-sectional area of power piston 
flow gain of input servovalve 
flow gain of power piston 
total volume of input piston 
total volume of power piston 
valve gradient of input servovalve 
valve gradient of power piston 
supply pressure (power piston) 
supply pressure (input servovalve) 
bulk modulus 

total mass of power piston and 
propeller 

total mass of input piston and 
valve rod assembly 

total flow pressure coefficient of 
input piston and servovalve 

total flow pressure coefficient of 
power piston 



57,67 in 

1102.0 in 

1.54 in 3 / sec/PSI 
0.11 in 3 /sec/PSI 
749.71 in 3 

14326.0 in 3 

I. 571 in^/in 

II. 938 in“/in 

400.0 PSI 

600.0 PSI 

100000.0 PSI 
69.8757 lb-sec 3 /in 

2 

4.9 lb-sec /in 

3 

2057.08 in /sec/in 

33426.54 in 3 /sec/in 
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IV. SYSTEM STATE EQUATIONS 



The system equations were first written in a state variable format 
[6]. Thus, in matrix notation, the system state equations became 

X(t) = A X ( t) + B U(t) 



Here X(t) was the matrix of system state variables and U(t) was the 

matrix of system input. From the overall system equations (9-12), Y(t), 

(t) , X ( t) and P (t) were chosen as the basic problem variables, and 
36a m 

X^ and became the input variables. Thus, the applicable state 
variables were 



Y = X(l) 


/ \ 

X 

II 

X 


P m = X(3) 
m 


ha ' X<5 » 


Y = X (2) 


X = X(5) 


p m = X ( 3) 
m 


ha ‘ X<6 » 



Y = X (2) X = X (5) 



For example from equation (9) one can write 



X(2) = -p X (3) - -t- n 
M Jl M * 2 



(13) 



In a similar way, the system equations were rewritten in matrix 
notation (Figure 5) . The accompanying signal flow graph for the over- 
all system equations is presented in Appendix A. The state equations 

were solved for various servovalve stroke values, X , and for various 

va 

loading, F^, using the IBM simulation language CSMP III as shown in 
the next section. 
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Figure 5„ System Equations (Matrix Notation) . 



V. ANALYSIS OF SYSTEM PERFORMANCE 



Manufacturer's data was consulted in order to formulate a power 
piston velocity specification. The maximum power piston velocity 
required was a value of 0.5 in/sec. Common operational velocities 
ranged between 0.4 - 0.5 in/sec. 

A. PROCEDURE 

1 . Inputs 

Electro hydraulic servovalves have strokes ranging from 0.005 

- 0.010 in. for small valves up to 10 gpm capacity to 0.015 to 0.030 in. 

for 50 gpm valves. In this case, a 15 gpm valve was used and input, 

X , was set from 0.005 to 0.012 in. 
va 

An external loading, F^, was assumed to range from 0.0 to 
40,000. lbf in the simulations. This range of loading could be expected 
in operational conditions. 

2 . Outputs 

The scheduling of control parameters is shown in Figure 6. 

Y and X signified the displacement of the power piston and the input 
piston valve rod assembly. P^ was the load pressure at the power 
piston, and was the fluid spring force at the power piston. 
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Figure 6. Scheduling of Control Parameters. 
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B. COMPUTER SIMULATION RESULT 



As mentioned previously, computer simulation results were taken for 

valve stroke displacements, X , from 0.005 to 0.012 in. Figure 7 

va 

shows the resulting power piston displacement. As can be seen, the sys- 
tem response is good. 

Figure 8 portrays the velocity of the power piston. These data 
demonstrate that the velocity of the power piston is directly propor- 
tional to the displacement of the valve input, The natural 

hydraulic frequency of the system was approximately 95 Hz for the con- 
trol unit input servovalve and piston and 110 Hz for the power piston 
assembly. These calculations are shown in Appendix B. 

Figure 9 is a plot of the load pressure, P , for various external 

loads. Without external loading, P was almost ze,ro . After an added 

m 

external load, P was changed slightly for a range of F n values from 

m 36 

2000.0 to 40,000.0 lbf. 

Calculations show that the displacement of the power piston upon 

# 

loading is 0.0016 in. for 40,000.0 lbf. Compared to the maximum move- 
ment (10.0 in.), this difference is about 0.016 percent. Thus, the 
displacement of the power piston was almost independent of external 
loading . 

The bulk modulus is a fluid parameter which characterizes the 
compressibility of the liquid. As can be seen from Figure 10, the 
value of the bulk modulus has a slight effect on power piston response 
time. Between bulk modulus values of 100,000 and 300,000 PSI , the 
power piston response was changed from 0.56 to 0.34 sec. Compared to 
a maximum required time of 19 to 20 sec, during which the power piston 
moved 10.0 in., this difference in power piston response time is very 
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Figure 7. Displacement Response of Power Piston 
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Figure 8. Steady State Velocity Response of Power Piston. 
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Figure 9. Pressure Response of Power Piston. 




Figure 10. Bulk Modulus Effects on Response. 
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small. This is an encouraging result since minimizing entrapped gas to 
insure high 6^ values is a difficult and costly process. 

Figure 11 shows the influence of bulk modulus on power piston 
velocity response. Again, in view of a total system response time 
requirement of 19 to 20 sec, the initial oscillations that occur within 
the first 0.3 to 0.6 sec of system operation can be considered negligible. 

Finally, in Appendix C is shown the computer program developed in 
this study and a typical computer output (Figure 15) . 



32 



300,000. PSI 
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Figure 11. Initial Velocity Response of Power Piston 



VI. CONCLUSIONS 



In this study, several types of C.R.P. systems were analyzed. These 
included: 

1. THE TWO VALVE -PISTON COMBINATION TYPE WITH SOLID CONNECTING ROD 
Results for this system have been presented in the previous 
sections of this report. 

2. A TWO VALVE-PISTON COMBINATION TYPE WITH CONNECTING FLUID LINE OF 

a. circular pipe construction 

b. annular pipe construction 

c. double annular pipe construction 

In all of the studies involving fluid lines, the fluid lines had very 
high frequency and very low damping ratio characteristics. Thus, essen- 
tially the influence of the connecting lines on response could be neglec- 
ted. Appendix D shows the result of various pipe line response calcula- 
tions . 

Computer simulation proved to be a useful design tool regarding the 
C.R.P. system. Full scale prototype testing, however, is the only con- 
clusive method of proving C.R.P. system response. Unfortunately, testing 
is both costly and time consuming. Thus, the computer simulation 
developed in this study could be used to designate critical testing 
instances and to identify those areas where design efforts would be most 
productive. In this way, the time involved in optimizing system perform- 
ance as well as subsequent production costs could be reduced. Likewise, 
the present computer simulation could be useful in specifying the correct 
approach for design testing practices. 
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APPENDIX A: SYSTEM SIGNAL FLOW GRAPH 
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Figure 12. System Signal Flow Graph. 



APPENDIX B: SYSTEM RESPONSE 



The governing equations could be separated into two parts - 

part (servovalve andpiston) and the output part (power piston) . 

Appendix, the hydraulic damping coefficient, £, and the natural 

to , of each part are determined, 
h 



A. INPUT PART (SERVOVALVE AND PISTON) 

From the Figure 13, the characteristic equation is 



C.E = S 



46 K 46 A 

V } MV 



2 2 
= S + 2 £ (0 S + (0, 
h h 
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Figure 13. Signal Flow Graph of Input Parts. 
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For nominal values of bulk modulus, the following results 



pertain: 



6 (psi) 

e 


co h (Hz) 


c 


100,000 


95.77 


0.00488 


300,000 


159.83 


0.00845 



B. OUTPUT PART (POWER PISTON) 

From the Figure 14, the characteristic equation is 



C.E = 



, 48 K _ 48 A 

s 3 + c — a-as.) s 2 + (— 2-E-) S 



V 



m 



V 

£ m 



48 K A 
e qm m 

V M/ 

m £ 



2 2 
= (S + R) (S + 2£ CO S + co, ) 

h h 



Again, for nominal values of bulk modulus, the following results 
were determined: 



8 (PSI) co, (Hz) 

e n 

100,000 110.45 

' 300,000 191.29 



c 

0.00897 

0.04094 




Figure 14. Signal Flow Graph of Output Parts. 
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APPENDIX C: COMPUTER PROGRAM 



INITIAL 

CONSTANT BE= 100000.0 
CONSTANT FF=2000 . 0 
CONSTANT CEM=1.54 
CONSTANT QMK=33426.54 
CONSTANT QAK=2057.08 
CONSTANT XM=69.8?57 
CONSTANT VM=4.9 
CONSTANT VA=749-71 
CONSTANT VC=14326.0 
CONST A ITT AM= 1102.0 
CONSTANT AA=57.67 
CONSTANT W=1 1.938 
CONSTANT PS =400.0 
CONSTANT PSA=600 . 0 

CONSTANT T1=0.0,T2=0.0,T3=0.0,T4=0.0,T5=0.0,T6=0.0 
DYNAMIC 

PARAMETER Z=(0. 005, 0.0 1.0. 012 ,0.013) 

FC=FF*STEP(22.5) 

XI N1=Z*STEP( 0 .0) 

XI N2=Z*STEP( 5* 0) 

XIN3=Z*STEP( 10.0) 

XIN4=Z*3TEP( 12.5) 

XIN5=Z*STEP(15.0) 

XI N6=Z*3TEP( 20 . 0) 

XIN7=Z*STEP(25.0) 

XI N 8 =Z*STEP( 30*0) 

XVA=XI Nl-XI N2-XI N3+XI N4-XI N'5+XI N6+XI N7-XI N8 

ydd=(am/xm)*pm-(i./xm)*fc 

yd=intgrl(ti,ydd) 

Y=INTGRL(T2,YD) 

pmd=-(4.*be*qmk/vc)*y-(4.*be*am/vc)*yd-(4.*be*cem/vc)*pm . . . 

+(4.*be*qmk/vc)*x 

PM=INTGRL(T3,PMD) 

fk=o.86*w*(ps-pm)*(x-y) 

xdd=(fk/vm)*y-(fk/vm)*x+(aa/vm)*pla 

xd=intgrl(t4,xdd) 

X=INTGRL(T5,XD) 

pud=-(4.*be*aa/va)*xd-(4.*be*cea/va)*pla+(4.*be*qak/va)*xva 

pla=intgrl(t6 5 plad) 

OUTPUT Y 
PAGE MERGE 
METHOD RKSFX 

TI MER FI NTI M=35 . 0 , 0UTDEL=0 . 5 , PRDEL=0 . 5 , DELT=0 . 002 

END 

STOP 
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Figure 15. Computer Program Output 



APPENDIX D: FLUID LINE RESPONSE 



When a fluid is used for transmitting pressure variations, it is 
necessary to determine the response of the transmission system. In 
the course of this study, the following line response determinations 
were formulated • 
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Figure 16, Circular Type of Line Response Determinations. 
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Figure 17. Annular Type of Line Response Determinations. 
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Figure 18. Double Annular Type of Line Response Determinations. 
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